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ABSTRACT 



O 

CD , 

Q I We review how to obtain the thermodynamic Bethe Ansatz (TBA) 

equations for the antiferromagnetic Heisenberg ring in an external 
^ ■ magnetic field. We review how to solve these equations for low tem- 

perature and small field, and calculate the specific heat and magnetic 
^ ! susceptibility. 

(N 
(N 

1. Motivation 

(N 

a^ 

^ I Quantum spin chains provide some of the simplest examples of physical systems 

with quantum-algebra^ symmetry. A well-known example is the open quantum spin 
D I chain consisting of N spins with the Hamiltonian 



N-l 



k=l 



k"k+l 



l{Q-Q-'){^l~cr%), (1.1) 



where a are the usual Pauli matrices, and q is an arbitrary complex parameter. 
This model is'^'"^ integrable and the Hamiltonian commutes^'^ with the generators 
of the quantum algebra Uq[su{2)]. 

This model can be generalized in various directions. We have considered in- 
tegrable open chains with quantum-algebra symmetry that have instead spins in 
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higher-dimensional representations , or spins in the fundamental representation of 
larger symmetry algebras^. (See also Ref. 8.) 

We would like to determine the physical properties (e.g., low-temperature mag- 
netic susceptibility and specific heat) of such models for |g| = 1 in the thermody- 
namic {N — > cx)) limit. For the model (1.1), some steps in this direction have 
already been taken in Ref. 4. Because these models are integrable, the spectra 
of their Hamiltonians are known implicitly through the Bethe Ansatz (BA) equa- 
tions. Unfortunately, this is still quite far from determining the models' physical 
properties. 

We consider here instead, as a warm-up exercise, the analogous problem for 
the antiferromagnetic Heisenberg ring, with Hamiltonian 



where H{> 0) is an external magnetic field. Although this model is comparatively 
simpler, the methods used to investigate it can presumably also be used to study 
integrable chains with quantum-algebra symmetry. 

The low-temperature thermodynamics of this model was determined approxi- 
mately twenty years ago as the result of the work of many people: Yang and Yang^, 
Takahashi^'^, Gaudin^^, Johnson and McCoy^^, et al. The fact that this work (cer- 
tain details of which have never been published) is scattered among many papers 
contributed to the difficulty we experienced in learning about it. 

We summarize this work here, in the hope of making it more accessible to 
others. Following Takahashi, we formulate in Section 2 the so-called thermody- 
namic Bethe Ansatz (TBA) equations for the antiferromagnetic Heisenberg ring. 
These are an infinite set of coupled nonlinear integral equations, which describe 
the equilibrium thermodynamics of the model at temperature T and in the field 
H. Following Johnson and McCoy, in Section 3 we solve the TBA equations in a 
systematic low-T and small- perturbative expansion; and we calculate the free 
energy, specific heat, and magnetic susceptibility. 

We warn the reader that we shall encounter along the way several arguments 
that are far from rigorous. Nevertheless, as far as we know, this is still the present 
state of the art. 

While there are other ways^^'^^ of calculating the low-temperature specific heat 
within the general TBA approach, the method pursued here has the merit of treating 
this calculation in the same manner as the one for the magnetic susceptibility. 





(1.2) 
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Review articles on the thermodynamics of integrable models which we have found 
particularly useful are Refs. 15 - 18. 



2. TBA equations 

The Hamiltonian for the antiferromagnetic Heisenberg ring in a magnetic field is 
given by Eq. (1.2). The corresponding energy eigenvalues, which can be determined 
by either the coordinated^, algebraic^°, or analytic^ -"^ Bethe Ansatz, are given by 

where the variables {Aq,} satisfy the Bethe Ansatz (BA) equations 
These equations can be written more compactly as 

M 

[ei(A,)]^ = -ne2(Aa-A^), (2.3) 

13=1 

where 

\ I in 

^n(A) ^ . (2.4) 

These equations are difficult to solve for finite N. Fortunately, in the physically- 
interesting thermodynamic (A^ oo) limit, these equations become easier to solve. 
Indeed, one then finds (see, e.g., Ref. 15) that the BA equations have "string" 
solutions. For instance, a string of length 1 is 

AW=a;, 

where x is an arbitrary real number, which is called the "center" of the string. 
Moreover, a string of length 2 is 

A(i) = x + I 
A(2)=a;-| ' 



a string of length 3 is 




and so on. We adopt the "string hypothesis" which states that all the solutions 
{Aq , a = 1, • ■ ■ , M} organize themselves into such strings. That is, the solutions 
are collections of M„ strings of length n of the form 

XinJ) = Xl + ^(!^-j) , (2.5) 



where j = 1, ■ ■ ■ , n; a = 0, 1, ■ ■ ■ , M^; n = 1, ■ ■ ■ , oo; and the centers AJJ are real. A 
particular solution of the BA equations corresponds to a set of non-negative integers 
{M„} and the M„ real numbers AJJ for each n. Observe that the total number of A 
variables, and BA equations, is M = nM^. 

Under this hypothesis, the BA equations (2.3) become 

oo Mm m 

ei(A(r'^))] = - n n n ^^(a^^'^ - >^^r'^) • (2.6) 

m=l/3=0fc=l 

Following Takahashi-^° and Gaudin^^, we form the product nj=i; obtain a set 
of equations for the (real) centers A^: 

00 Mm 

[en{K)f = i-r n n ^nmiK - A^) , (2.7) 

m.=l ^=0 

where 

-Enm(A) = e|„_^|(A) ef„_^|+2(A) • • • e^+^_2(A) en+m(A) . (2.8) 
In obtaining this result, we use the lemmas 

n min{n,rn) 

n^-(^^'^)= n e„+^+i-2KA2), (2.9) 

i=i 1=1 

n m 

n n ^2(Ai"'^'^ - A^"^'')) = £;„^(A2 - A^) . (2.10) 

j=lfc=l 
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Since the equations (2.7) involve only products of phases, it is useful to take 
the logarithm: 

oo Mm 

Npr^iK) = 27, + E E -nrniK - A^) , (2.11) 
m=l/3=0 

where 

Pn(A) = ^lnen(A) , S„^(A) = zln£;„^(A) , (2.12) 

and are integers or half-integers constrained by 

~'^max — — '^max ■ (2.13) 

Faddeev and Takhtajan have given-*^^ a prescription for determining an explicit 
expression for J^^j. from (2.11). Any integer or half-integer in the range (2.13) is 
called "admissible". We adopt the important hypothesis that the numbers {J^} 
can be regarded as quantum numbers of the model: for every admissible set {J^} 
(no two of which are identical), there is a unique solution {A^} (no two of which 
are identical) of Eq. (2.11). Such solutions are called "particle rapidities". 

As we shall see, it is useful to introduce (following Yang and Yang^) the concept 
of "hole" . Given a set of quantum numbers {,7^}, we define the set of numbers { J^} 
to be the admissible values not in { J^}. We define corresponding quantities {A^}, 
called "hole rapidities" , by 

oo Mm 

Npr^iK) = 27rJ: + E E ^-miK - A^) . (2.14) 

m=l /3=0 

For N ^ oo, the distributions of the numbers {A^} and {A^} along the real 
axis are characterized by corresponding particle and hole densities Pn(A) and p^(A), 
respectively: 

Npn{X)dX = number of A'^'s in dX 

Npn{X)dX = number of A'^'s in dX . (2.15) 

We must now reformulate the description of the model in terms of these densities. 

We begin by exhibiting a constraint between the hole and particle densities. 
To this end, we define the function 

( oo Mm 1 

h-{X) ^^{Npn{X)-J2Yl ^r.m{X ' Xf) \ . (2.16) 
I m=l (3=0 I 



5 



'a 5 



Evidently, 

h-CK)=Ja- (2.17) 

We make the assumption, which is imphcit in Refs. 10 and 11, that h'^{X) is a 
monotonic increasing function of A. It follows that 

N [p„(A) + p„(A)] dX = number of A'^'s and A'^'s in dX 

= number of times /i"(A) ranges over J^'s and J^'s in dX 
= /i"(A + dA) -/i^(A) = d/i^. (2.18) 

That is, 

1 dh"^ 

Pn{X)+Pn{X)^-^. (2.19) 

Making in Eq. (2.16) the following replacement of sums by integrals 

^E()^/ {)PmiX')dX', (2.20) 

we see from (2.19) that 

m=l ^ 

This can be rewritten in the final form 

oo 

Pn + ^ ] ^nm * Pm — 0,n j (2.22) 

m=l 



where 



and 



= ^{X)5nm + (1 - ^nm) «|n-m| W + a„+^(A) 
min{n,m) — l 
+ 2 ^ a| 

n — m[+2l 

(A) , (2.23) 



,^(A)^J-%i^ = J-^L^, (2.24) 
^ 27r dA 27r A2 + if ^ ^ 
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and * denotes the convolution 

/oo 
dX' fix - X')9{X') = ig*f) (A) . (2.25) 
-oo 

One can regard the constraint equation (2.22) as a definition of pn in terms of pn- 
We introduce here for future reference the "inverse" function 

A-^{X) = 5{X)5nm - S{X) {5n,m+l + 5n,m-l) , (2.26) 

where 6(A) is defined by 

s{X) = — ^ . (2.27) 
^ ^ 2ch7rA ^ ^ 

It has the properties that 

oo 

J2 iKn' * ^n'm) (A) = S{X)6nm , (2.28) 
n' = l 

and 

oo oo 

E (^nm * «m) (A) = s{X)6nl , J] * m = . (2.29) 

m=l m=l 

These relations can be easily derived with the help of Fourier transforms, for which 
we use the following conventions: 

/oo -1 pco 

e-^ f{X)dX , /(A) = — y_ e--^ f{aj)doj . (2.30) 

In particular, we note the convolution theorem 

/oo 

e'^^ if * g) (A) dX = f{u) g{u) , (2.31) 
-oo 

and the identities 

«n(a;) = Q ^ ^ , .(c) = ^ , (2.32) 
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An^(a;) = (^cth ^ j ^e-|n-m||a;|/2 _ g-(n+m)|a;|/2j ^ ^g.SS) 

^nm('^) = ^nm - K^) (^n.m+l + 5n,m-l) • (2.34) 
7 



We recall that the Bethe Ansatz provides the expression (2.1) for the energy- 
eigenvalues in terms of A's. In order to recast this expression in terms of densities, 
we first invoke the string hypothesis to arrive at 



^=-^EEE7— ^^^1 H[^-j:nM„]. (2.35) 



We next observe that 



n ^ . , n . J n 



4 

= ^dXs ^""'"^^^^ = ^^Pn(AS) =7ra„(AS), (2.36) 

where, in passing to the second line, we have used the lemma (2.9). It follows that 
the energy is given by 

oo M„ / oo 

n=l a=0 \ n=l 



= -iVTT 



oo /"OO / -1 oo /"OO \ 

E / On(A)p,(A)dA -NH Pn{X)dX . (2.37) 

n=l-^-~ \^ n=l -^-"^ / 



Given a set of particle densities {pn}, this expression determines an energy level of 
the system. Notice that this expression does not depend on the hole densities pn- 

As an illustration of this formalism, let us briefly consider the ground state of 
the system in zero field {H — 0). We assume that for this state, there are no holes. 
That is, 

Pn = . (2.38) 

(We shall prove this later from the TBA equations.) It follows from the constraint 
(2.22) that 

oo 

E "^nm * Pm = an- (2.39) 

m=l 

Thus, the particle densities are given by 

■oc 

Pn = E ^"'^ * am = s6nl , (2.40) 
m=l 
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where s(A) is given by (2.27). That is, the antiferromagnetic ground state is de- 
scribed by strings of length 1 ( i.e., real solutions of the BA equations). From 
Eq. (2.37), we conclude that the energy per site for the ground state is 



eo 



El roo 

-^ = -7r/ ai(A)s(A)cZA = -ln2, (2.41) 



which is the well-known result of Hulthen^^. 

This example raises two important questions: 

(i) how to establish the result (2.38) that the ground state has no holes? 

(ii) how to study low-lying excitations? 

An answer to both questions is to consider the system at temperature T > 0. The 
equilibrium state at T = is the ground state; and the equilibrium state at low but 
nonzero temperature provides information about the low-lying excitations. 

We now argue heuristically, following Refs. 9, 17 and 18, that the equilibrium 
state at temperature T is described by certain equilibrium densities p^^ and p^. 
We begin by observing that the partition function is given by 

Z = tr e-^/^ = ^'"^'^ = E ^ ^Xip, P)] W[p, p] e-^H/^ , (2.42) 

eigenstates p,p 

where the delta function enforces the constraint 

X{P,P) = 0, (2.43) 

which represents the constraint (2.22), and W[p, p] is equal to the number of states 
corresponding to given densities p and p. Introducing the entropy 

S = \nW, (2.44) 

and using the fact that the free energy is given by 

F = E-TS, (2.45) 

we see that 

Z = J2^ [X(P, P)] e-(^[^l-^^[^'^l)/^ = J2 ^ [X(P, P)] e-^[^'^]/^ . (2.46) 
p,p p,p 

Using the saddle-point approximation (for M , N ^ oo with M/N — constant), we 
conclude that 

Z = e-^^P-^'P-^y^ , (2.47) 
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where the equihbrium densities are determined by the condition 

5F = , (2.48) 

P=Peq ,P=Peq 

subject to the constraint (2.43). 

In order to compute the entropy (2.44), we observe (again, foUowing Yang and 
Yang^) that 

dSn = In {number of states (strings of length n) in dX} 
= In {number of ways of choosing A'^'s in dX} 



= In 



[N{pn + Pn)dX]\ 



{NpndXy. {NpndXy. 

~ N{{pn + Pn) In (pn + Pn) " Pn In 

Pn Pn 

\npn}dX, (2.49) 

where, in passing to the last line, we have used Stirling's approximation. We con- 
clude that the entropy is given by 

^ = ^ dSn^N^ {{pn + Pn)^n{Pn + Pn) - Pn^nPn- Pn^nPn}dX. 

n=l'^-°° n=l"^-°° 

(2.50) 

The hole concept was introduced precisely for calculating the entropy. 

The free energy F of the Heisenberg ring is given by (2.45), with the energy E 
given by (2.37) and the entropy 5" given by (2.50). The condition dF = subject 
to the constraint (2.22) yields the thermodynamic Bethe Ansatz (TBA) equations 

oo 

Tin (^1 + e^"/^) ^J2^nm*Tlii(l + e"^-/^) - Tron + nH , (2.51) 



m=l 



where the functions e„(A) are defined by 



ptm 

Pn%X) 



6n(A)=Tln(^) . (2.52) 



The TBA equations are an infinite set of coupled, nonlinear integral equations, 
which we shall assume to have a unique solution {e^(A)}. It follows that {en{X)} 
are even functions of A, 

6n(-A) = e,(A). (2.53) 

Having solved for these functions, one can then use the constraint equations (2.22) 
to determine the equilibrium densities {p^'^(A)} and {p^'^(A)}. 

10 



The free energy per site at equilibrium is given by 

^ = -TY.j JX an(A) In (l + e^-^^)/^) - -H . (2.54) 



n=l ' 



This result is obtained by using the constraint equations (2.22) to eliminate pn from 
the expression for F/N, and then noting that the coefficient of p„ vanishes by virtue 
of the TBA equations. Alternatively, by eliminating p„, one can show that 



dX s(A)ln 1 + , (2.55) 

-00 



where cq is given by (2.41). This expression for the equilibrium free energy is 
particularly useful, since it depends on only one of the en(A)'s. 



3. Thermodynamics for small T and H 

In this section, we shall solve the TBA equations (2.51) for ei(A) for small 
values of T and H, and evaluate the expression (2.55) for the free energy. 

As we have seen, the TBA equations are an infinite set of coupled equations. A 
significant simplification occurs for T — > 0: namely, the equation for ei(A) becomes 
decoupled from the other equations. In order to see this, we first cast the TBA 
equations in the form given by Takahashi^°: 

ei = -Trai +H + a2*Tln(l + e"'^/'^) 

+ (ao + a2)*^am*Tln(l + e-^-+i/^) , (3.1) 

m=l 

en = H + ai*T\n(l + e'^-^/^) + 03 * Tin (^1 + e"'"/^) 

00 

+ (ao + a2)* ^a^*Tln(l + e-''"+"/^) , n>2, (3.2) 

m=l 

where ao(A) = S{X). To obtain Eq. (3.1), we start by writing the n = 1 TBA 
equation using the explicit form of as given in (2.23). We then rearrange the 
terms and use the identities 

To obtain the equation for €2, we use the additional identity 

an*H = H (3.4) 
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to write the driving term —7ra2+2H of the n = 2 TBA equation as ai*(— Trai + H) + 
H, and we then use Eq. (3.1) to ehminate the terms in parenthesis. The other 
equations are obtained by following the same strategy. 
From this form of the equations, we conclude that 

en(A)>0 forn>2. (3.5) 

This positivity property will play an important role in the arguments that follow.* 

Consider now Eq. (3.1). The positivity property (3.5) implies that exp(— e^+i/T) 
goes to zero exponentially as T — > for m > 1 and can therefore be neglected. This 
is the reason for the decoupling of equations mentioned above. On the other hand, 
ei(A) can have either sign. Defining e{X) to be the T — > limit of ei(A), we see that 

lim Tin fl + e^'^^^) = ±£^ , 

where we use the standard notation 

S~ = ^{S- \e\) , £+=£-£ 

It follows from (3.1) that 

£ = — TTOi + H — a2* £~ ■ (3.8) 

It will be convenient for subsequent analysis to work with the following equiv- 
alent equation, which does not involve e~: 

e = — TTS + /i * £+ , 

where s(A) is defined in (2.27), and h{X) is defined by 

/i = s * ai . (3.10) 

In order to obtain (3.9), we substitute e~ = £ — £+ into (3.8); and we then solve for 
£ in terms of £+ with the help of Fourier transforms. We observe for future reference 

* An alternative way of demonstrating this positivity property relies on the "in- 
verted" TBA equations, which are obtained by forming the convolution of the TBA 
equations with A~}^. However, these equations must be supplemented with a con- 
dition on €n for n — > cxD. 
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(3.6) 



(3.7) 



(3.9) 



that by a similar procedure, (3.1) can be recast (after neglecting the infinite sum, 
but before taking the T — > limit of ei) as follows: 

ei = ^H -7TS + h*Tlii(l + e'l/^) . (3.11) 

Evidently, the T ^ limit of this equation gives (3.9). 

We have obtained the linear integral equation (3.9) for e. Of central importance 
in obtaining this result was the positivity property (3.5). We shall now show that 
Eq. (3.9) can be solved, and that it has a unique solution. We shall also establish 
certain properties of the solution, which will turn out to be crucial for formulating 
the linear integral equation obeyed by the leading correction to this lowest-order 
result. 

We first briefly consider the case H = 0, which corresponds to the ground 
(vacuum) state. Since s(A) is positive for all A, Eq. (3.9) has the solution 

.(A) = -..(A) = -5^. (3.12) 

From the definition (2.52), the constraint equation (2.22), and the positivity prop- 
erty (3.5), it follows that the ground-state densities of holes and particles are pre- 
cisely those anticipated in the previous section. (See Eqs. (2.38) and (2.40), respec- 
tively.) 

We turn now to the case H ^ 0. We observe, following Yang and Yang and 
Takahashi, that Eq. (3.9) can be solved by iteration. Indeed, set 

£(0)(A) = ^ff-7rs(A), (3.13) 
£(n)(A) ^^H- ns{X) + (h * (A) , n > 1 . (3.14) 

If the sequence {£(o)(-^) : r "} has a limit 

lim £(„)(A) =£(A), (3.15) 

n— >cxD 

then this limit satisfies Eq. (3.9). In order to show that this sequence is convergent, 
we prove that it is monotonically increasing (i.e., e(„)(A) > (A)) and that it is 

bounded from above. To demonstrate that the sequence is monotonically increasing, 
we proceed by induction, using the fact that h{X) is positive for all A, and that 

£(n) = £(n-l) + h * - £+ _2)) , n>2. (3.16) 
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To show that the sequence is bounded from above, we cast the relation (3.14) in 
the form 

£(n) = £(0) + ^ * £(0) H h /i * • • • * * £+p (3.17) 

where in the last term the factor h appears n times. From (3.13), we see that 
£(o)(A) < H/2. It follows that 



H 



and therefore 



We conclude that 



h*---*h*e+^{X) < ( - ) — , (3.18) 



on+l _ 1 

^in)W<^^H. (3.19) 



£(n)(A)<iy (3.20) 



for any n, and so the sequence is bounded from above. This completes the proof 
that the iteration procedure (3.13), (3.14) gives a solution of Eq. (3.9). 

From this iterative solution, we can understand the general behavior of £(A) 
as a function of A. For A > 0, this function is monotonically increasing. Indeed, 
following again Refs. 9 and 10, we prove the monotonicity of £(„)(A) by induction. 
From (3.13), we determine that £'(o)('^) > for A > 0. Differentiating (3.14) with 
respect to A, we obtain 

/OO 
- ^tn-i)(0 dC (3.21) 
-OO 
/•OO 

= <o) (A) + [h{X-0-h{X + C)] £+'_!) (0 dC . (3.22) 

Since the function h{X) is even and monotonically decreasing for A > 0, the expres- 
sion inside the brackets in Eq. (3.22) is positive. Therefore, if £(jj_i)(A) > for 
A > 0, then £(^-)(A) > for A > 0. This concludes the proof that £(„)(A) (and hence 
£(A)) is a monotonic increasing function of A, for A > 0. 

On the other hand, we have seen that £(A) is bounded from above. Therefore, 
the function £(A) has the following general behavior: 

figure available on request 

In particular, £(A) has only one zero for A > 0, 

£(a) = 0. (3.23) 
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Hence, s{X) < for A in the interval (0,0;), and £(A) > for A in the interval 
(a , oo). 

Observe that Eq. (3.9) is temperature-independent. Using the same approxi- 
mations in the expression (2.55) for the free energy this too will be temperature- 
independent. As emphasized by Johnson and McCoy^^, to find the leading order 
temperature dependence one must compute the leading correction to the solution 
ei = £ of the linearized equation (3.9). In order to obtain this correction, we make 
the substitution 

ei(A)=£(A) + 77(A) (3.24) 

in (3.11) and expand to leading order in rj. Since e is a solution of the linearized 
equation (3.9), we shall find an inhomogeneous term in the resulting equation for 
rj. Indeed, we have that 



77 = * jrin [1 + e(^+^)/^] - £+} 



(3.25) 



Because £(A) is an even function of A with a single zero, at A = a, for positive A, 
we see that 



+ 



—00 

a 



dX 



' /i(A- A'){ 



Tin 



1 + e 



(e(A')+r,(A'))/T 



-S{X')] 



+ f dX' HX - A') Tin [1 + e(-(^')+r;(A'))/Tl 

J —a 

/—ct po 
-00 J a 



dX' h{X-X') ri{X')+Eh{X) 
where the inhomogeneous term Eh{X) is given by 



(3.26) 



(3.27) 



For T — > 0, the major contribution to the integral comes from the regions near the 
zeros of e, so we expand e(A) about A = a. 



s{X) = t{X-a)^0{{X-af) , 

One then finds that the leading T- dependence of is 

2T2 



t = 



de 
dX 



(3.28) 



Eh{X) = ^ [h{X -a) + h{X + a)] / du In (l + e"") 

Jo 

[h{X -a) + h{X + a)] . 



t 



6t 



(3.29) 
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Thus, we obtain the following linear, but nonlocal, equation for 77,* 



r;(A) =U +J \d\' h{\- A') r]{\') + [h{\ - a) + h{\ + a)] . (3.30) 

This equation for r] and that of (3.9) for e complete our results for the T - expansion 
of ei. We turn now to the expansion for small values of by which from Eqs. 
(3.9) and (3.30) we shall generate a pair of integral equations of the Wiener-Hopf 
type. (For a detailed exposition of Wiener-Hopf equations, see Ref. 23.) 

For iy — > 0, we see from the solution of the zeroth-order equation (3.13) that 
the function e{\) has a zero for A = a = 0{\n.H). Assuming an expansion in 
powers of (Inff)"^, we conclude that 



1 

a = 

TT 



In ( —H] + \nK + 0' ^ 



V27r y \\nH 



(3.31) 



where the constant k (which is independent of H) has still to be determined. 
It will prove convenient to work with the functions 

5(A) = {^'"''f + A>0 _ (3 3^j 

and 

T(A) = I StS^(^ + ^ > , (3.33) 

' \ A < ' ^ ' 

instead of the functions £(A) and ri{X). The factors e""' and e~^"/T^ in (3.32) and 
(3.33), respectively, are chosen such that the driving terms in the equations for S{X) 
and T(A) have a nonvanishing limit as T and i7 — » 0. 

We return now to Eq. (3.9). Following Refs. 24 and 12, we write the limits 
of integration explicitly, keeping in mind that e{\) > for —00 < A < —a and for 
ct < A < 00; and we shift the integration variables so that they run from to cxd. 
Since h{\ + 2a) vanishes as — > (ct — > 00) for finite A, we see that the i7 — > 
limit of the equation for S (A) is 

POO 

S{X) ^tt{1- Ke-""^) + dX' h{X- X') S{X'), A>0. (3.34) 

^0 



* While Johnson and McCoy obtain this equation from the "inverted" TEA equa- 
tions, we obtain it directly from the TEA equations. 
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Similarly, from (3.30), we see that the — > limit of the equation for T(A) is 

T(A) = + dX' hiX - A') T(A') , A > , (3.35) 

where 

S'{0)=^S{X) ^e'^'^Kt. (3.36) 



dX 



A=0+ 



These equations can be written in the standard Wiener- Hopf form 

/oo 
dX' h{X - X') S{X') , 
-oo 



-oo 
foo 



where 



/oo 
dX' h{X - X') T(A') , 
-oo 

— oo < A < oo , (3.37) 



f (\\-) 7r(l-«e-'^^) A>0 . . 

/s(A) - <! Q A < ' ^^-^^^ 



bs{X) - ^ _ J-^ /i(A - A') 5(A') A < ' ^^"^^^ 



and similarly 



MX) = I S'(0) ^ > U , (3.40) 

6t(A) = I _ ^oo^ ^ _ y^^,^ ^ > Q . (3.41) 

We shall solve these equations by Fourier transform. Since S{X) and T(A) 
vanish for A < 0, the Fourier transforms S{uj) and T{lj) are analytic in the upper- 
half- plane Im a; > 0, which we denote by 11+ . Observe that 



2 



oo 



S{Q) = — duj S{u) , (3.42) 

where the factor 2 in the numerator is due to the discontinuity of S{X) at A = 0. 
Evidently, the contour integral of <S'(a;) along a closed contour in 11+ is zero. By 
explicitly evaluating the contribution of the semi-circle at infinity, we obtain 

I duj S{oj) = - I duj S{u}) = -i7T lim ujS{uj) , (3.43) 

J — oo J semi— circle |a;|— >oo 
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where the limit is taken in 11+ . On the other hand, S{0) = 0, as follows from 
s{a) = 0. Therefore, 

,5(0) = -i lim u;S{u;) = . (3.44) 



a; — >oo 



Similarly, 

S'{0) = ^S{X) =- lim u'Siu) . (3.45) 

In deriving this result, we must again include a factor 2 due to the discontinuity 
(here, of S'(X)) at A = 0. 

The Wiener-Hopf equations in Fourier space are 

Siu) = fsiio) + bsiu) + h{u) Siu) , (3.46) 
f{io) = friio) + briuj) + h{uj) f{uj) . (3.47) 

Since h{—uj) = h{ijj) and 1 — h{uj) ^ for — oo < a; < oo, there exists^^ a 
unique factorization 

{l-h{uj)^ ^ = G+{uj) G-{uj), -oo<a;<oo, (3.48) 

where G+(a;) and G_(a;) are functions that are analytic and nonvanishing in the 
half- planes 11+ and n_, respectively, and are normalized by the condition 

G+(oo) = ^-(oo) = 1. (3.49) 

Moreover, the fact h{—\) = h{X) implies that (for oj in n_) 

G^{u;)^G+{-u;). (3.50) 

As we shall see, explicit expressions for G+ and G- are not needed to compute the 
free energy to the order in which we work. A similar phenomenon occurs in the 
Wiener-Hopf calculations of Yang and Yang^'^. 

These properties will be used to solve the Wiener-Hopf equations (3.46) and 
(3.47), to which we now turn. Using the factorization (3.48), the equation (3.46) 
for S can be rewritten as 

G+^ S = G- (/5 + bs) . (3.51) 

Observe that the left hand side is analytic and bounded in 11+ , whereas G-bg is 
analytic and bounded in H_. The term G-fs has a decomposition as 

Gjs = P- (C-Zs) + P+ (gJs) , (3.52) 
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where P± (^G-fs^ is analytic in Il±. This decomposition is uniquely specified by 

the requirement P± (^G-fs^ — for a; — > cxd in II-i-. Taking the P-f projection of 
(3.51), we have that 

S^G+P+ (G-fs) . (3.53) 

From (3.38) we compute that 

fs{^) = ^) , (3.54) 

where one is to take e — > at the end. The decomposition (3.52) of G-fs is then 
found by subtracting the residues of /s, i.e., 

G-iu;)fs{u;) =nil^-l—{G-iu) - G_(-^e)) - (G_(cc;) - G_(-7r^)) 

+ TTz I -^—G_ (-ie) i-TTi) \ . (3.55) 



Hence, 



S{uj) = -^^G+{uj)G-{0) - -^^G+{uj)G-{-m) . (3.56) 
<jj + ie <jj + ni 



The boundary condition (3.44) implies that the parameter n is given by 

K = G-{-7Ti)-^ G-{0). (3.57) 
Using this result we conclude that S{u!) is given by 

S{u;) = 7Ti (^—. G+{u;)G_{0) . (3.58) 



From Eq. (3.45) and the fact that G+{u) ^ 1 as |a;| ^ oo in 11+ , we find that 

S'{0) = 7T^G-{0). (3.59) 

We turn now to the equation (3.47) for T(u;). Proceeding as before, we use the 
factorization (3.48) to arrive at the formal solution 

f = G+P+ (g_/t) (3.60) 
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(cf. Eq. (3.53) ). The explicit calculation of /t is difficult, but can be avoided by 
the following trick. Consider the functions /(A) defined as 

/(A) = |i|. (3.61) 

where A ranges over the entire real line. From Eq. (3.40), it is evident that 

/t(A) = /+(A), (3.62) 

where 

III Ill f.u.f- 

(3.63) 

The Fourier transform of /(A) is readily computed: 

From the factorization (3.48), it follows that 

/+ + /- = (1-GI^G;');57|^- (3-65) 

After multiplying both sides of this equation by G- , we see that the P-f projection 
of is given by 

P+ (g_/+) = a - , (3.66) 

where a is a constant. Requiring the right hand side to vanish for |a;| — > oo in 11+ 
determines this constant to be 

We conclude from (3.60) that T{uj) is given by 

T(a;) = (G+(a;)-l)-^, (3.68) 

where 5"(0) is given by (3.59). 

To summarize: we have made the expansion (3.24) of ei(A), and we have 
changed in Eqs. (3.32), (3.33) from the variables £(A), 77(A) to the variables 5'(A), 
T(A), respectively. Using Wiener- Hopf methods, we have determined in Eqs. (3.58), 
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(3.68) the corresponding Fourier transforms S{(jij), T{uj) in the hmits T — > and 
iy — > 0. These expressions involve G+(a;) and G_(a;), which appear in the factor- 
ization (3.48). We shall now use these results to calculate the free energy. 

Substituting the expansion (3.24) of ei(A) into the expression (2.55) for the free 
energy (keeping in mind the discussion immediately following (3.24)), we obtain 



N='' 3t ^ 



dX s{X) [e{X) + ri{X)] , 



(3.69) 



where the ground state energy per site eo is given by (2.41). 

We aim for a double expansion of F/N to quadratic order in both T and H. 
To this end, we rewrite (3.69) in terms of S{X) and T(A). Consider first the 
term. For i? — > 0, we can make the approximation 



s{a) = e~ 



(3.70) 



and hence (remembering (3.36)) this term becomes 

3S'{0) ' 

We next manipulate the integral in the last term as follows: 



(3.71) 



/ dX s{X) [e{X) + r]{X)] = dXs{X + a 

J a Jo 



+ T^nx) 



K 



6e" 



. (3.72) 



Taking the a — > oo limit of this expression, we arrive at the following expression for 
F/N, 



F H^K 
iV='°" 2^ 



(3.73) 



Writing the functions S{X) and T(A) as the Fourier transforms of S{ui) and T{ui), 
respectively, and then performing both the A and uj integrals, we obtain 



F 

iV = "° " 2^- 



-B, 



with 



A = KS{7ri) , 



B = K 



S'{0). 



(3.74) 



(3.75) 
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The parameter k is given in Eq. (3.57); and using the property (3.50), this 
parameter can be reexpressed as 

At = G+(0)G+(7ri)~^ (3.76) 

We consider now A. Evaluating S{'Ki) using the expression (3.58), we obtain 

^=^G+(0)G_(0). (3.77) 

RecaUing the factorization equation (3.48) and the fact that 

h{ijo) = s{uj) ai{ijo) , (3.78) 

we conclude that ^ 

A=]-(l-h{Q)) =1. (3.79) 



2 

There remains to compute B. Evaluating T(7rz) using the expression (3.68), 
we obtain (after a crucial cancelation) 

B=^. (3.80) 

Recalling the expression (3.59) for S"(0) and using again the property (3.50), we 
see that 

B=\. (3.81) 

TT 

The expression for the free energy per site is therefore 

As foreseen above, this result was obtained without using explicit expressions for 
G+ and G_. It follows that the magnetic susceptibility and specific heat, to lowest 
order, are given by 



X 



(F 



1 



For a critical chain, the free energy per site is given by^^ 



T. (3.83) 

H 3 ^ ' 
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where c is the central charge and Vs is the velocity of sound. For the antiferro- 
magnetic Heisenberg Hamiltonian (1.2), Vg = 7r/2 (see, e.g., Ref. 15), and hence 
c= 1. 
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